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Abstract 


Bats are flying mammals distributed worldwide known to host several types of Coronavirus (CoV). Since they were reported as 
the probable source of spillover of highly pathogenic CoV into the human population, investigating the circulation of this virus in 
bats around the world became of great importance. We analyzed samples from 103 bats from two distinct regions in Brazil. 
Coronavirus from the A/phacoronavirus genus was detected in 12 animals, 11 from Sao José do Rio Preto—SP region and 1 from 
Barreiras—BA region, resulting in a prevalence of 17.18% and 2.56% respectively. The virus was detected not only in intestines 
but also in lungs and liver. Phylogenetic analysis based on nsP12 genomic region suggests that the sequences group according to 
host family and sampling location. Studies on the circulation of these viruses in bats remain important to understand the ecology 


and evolutionary relationship of these pathogens. 
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Introduction 


Coronaviruses (CoVs) are single-stranded positive sense 
RNA viruses that belong to the Coronaviridae family. The 
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family comprises four genera: Alphacoronavirus (a-CoV), 
Betacoronavirus ((3-CoV), Gammacoronavirus (y-CoV), and 
Deltacoronavirus (6-CoV) [1]. These viruses infect several 
species of animals and can cause hepatic, enteric, neurologi- 
cal, and acute and chronic respiratory infection [2]. In humans, 
they generally cause mild upper respiratory tract disease with 
low mortality [3]. In 2003 and 2012 respectively, the emer- 
gence of highly pathogenic 3-CoV SARS-CoV (severe acute 
respiratory syndrome) and MERS-CoV (Middle East respira- 
tory syndrome) revealed that this group of viruses can also 
cause severe respiratory illness in humans [4, 5]. 

Bats are widespread flying mammals that are distributed 
worldwide. About 1150 bat species have been described and 
179 were reported in Brazil [6-8]. Numerous viruses have 
been reported to infect bats, and these animals are being stud- 
ied as potential reservoir for several zoonotic diseases, includ- 
ing Coronaviruses [3, 9, 10]. The majority of Bat Coronavirus 
(BtCoV) have been identified in bats from Asia, Africa, and 
Europe [11-19]; however, CoV has also been reported in bats 
from South American countries including Brazil, Trinidad, 
Costa Rica, Ecuador, Panama, and Mexico [20-25]. 

Although there is no report of CoV pathogenesis in bats, 
studies suggest that SARS-CoV and MERS-CoV have a zoo- 
notic origin in bats [26]. Therefore, studying the circulation of 
BtCoV in bats around the world is of great importance to 
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detect potential sources of spillover of pathogenic viruses into 
the human population. Here, we investigate the presence of 
coronavirus in bats from two distinct regions in Brazil. 


Methods 
Samples 


Bats were collected in Sao José do Rio Preto region located in 
southeastern Brazil at the state of São Paulo (SP) and in 
Barreiras region, located in the northeastern part of the country 
at the state of Bahia (BA). The first is a transition zone be- 
tween the Atlantic Forrest biome, composed by an ecosystem 
of semi-deciduous forests, and the Cerrado biome character- 
ized by savannahs. While the second region is situated within 
the Cerrado biome. Both regions are under tropical climates 
with dry winters (Tropical Aw according to Köppen-Geiger). 
The animals were collected in urban and peri-urban locations, 
using mist nets and euthanized according to the Guide for the 
Care and Use of Laboratory Animals of the Institute for 
Laboratory Animal Resources Commission on Life Sciences 
National Research Council [27]. Also, collections were made 
under environmental licenses from the Brazilian Institute of 
Environment and Natural Renewable Resources (IBAMA), 
SISBIO no.: 21707-1 (Process: 02027.001957/2006-02); 
46190-1 and the project was approved by the Ethics 
Committee on the Use of Animals of the Institute of 
Biosciences, Letters and Exact Sciences (CEUA/IBILCE— 
Protocol 135/2016). Sex was determined based on Vizotto 
and Tadei, Miranda et al., and Laval [28-30]. Intestines, liver, 
and lungs were removed from bats from Sao José do Rio 
Preto—SP and intestines and liver from bats collected in 
Barreiras—BA. All samples were stored at — 150 °C until 
further processing. 


RNA Extraction and cDNA Synthesis 


Samples were thawed, immersed in TRIzol (Thermo Fisher 
Scientific), and tissues were homogenized mechanically by a 
polytron homogenizer (Marconi, Piracicaba, SP BR). 
Extraction was carried out according to manufacturer’s in- 
structions. RNA pellet was resuspended in water treated with 
DEPC (Sigma Aldrich), quantified by NanoDrop 2000 spec- 
trophotometer (Thermo Fisher Scientific), and stored at — 
150 °C. The cDNA synthesis was performed using High 
Capacity cDNA Reverse Transcription Kit (Applied 
Biosystems). 


Coronavirus Detection and Sequencing 


The presence of Coronavirus in the organs was tested by a 
pan-Coronavirus NESTED-PCR as described by Chu et al. 
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targeting the RNA-dependent RNA polymerase (RdRp), 
nsP12 region [31]. Positive samples were purified using 
DNA Clean & Concentrator (Zymo Research) and sequenced 
by Sanger method [32], using the set of primers from the 
NESTED reaction, in a 3130 XL Genetic Analyzer (Applied 
Biosystems). Sequences were submitted to Electropherogram 
Quality Analysis tool from Embrapa Genetic Resources and 
Biotechnology website (http://asparagin.cenargen.embrapa.br/ 
phph/). The platform performed the quality analysis using 
Phred [33] and contig assembly by CAP3 [34]. 


Sequence ID 


Sequences are identified by BtCoV, followed by the region 
where the host was collected (BA—Barreiras; SJRP—Saéo 
José do Rio Preto), the reference of the animal (e.g., 03, 04 
...) and by the identification of the organ where the viral RNA 
was detected (I—intestines; Li—liver; Lu—lungs). 
Nucleotide sequences were submitted to GenBank under the 
accession numbers MH974764 to MH974780, and specific 
identification is available in Online Resource 1. 


Datasets and Phylogenetic Analysis 


Two datasets were assembled to perform the phylogenetic 
analysis. The first one (dataset 1) was composed by all the 
sequences generated in this study plus sequences representing 
all four Coronavirus genera (Online Resource 1): a-CoV, 6- 
CoV, 7-CoV and 6-CoV. The second dataset (Dataset 2) is 
composed by our sequences and sequences available in 
GenBank of the nsP12 region of BtCoV from the 
Alphacoronavirus genus (Online Resource 2). Sequences 
were aligned using Muscle [35, 36], nested in the Seaview 
4.6.1 package [37]. The datasets were edited using BioEdit 
7.2.5 package [38]. 

A Bayesian Markov Chain Monte Carlo (MCMC) estima- 
tion of the phylogenetic trees based on the sequences in both 
datasets was done using BEAST v1.8 [39]. The analysis was 
carried out using General Time Reversible model plus gamma 
distribution (GTR + G) for dataset 1 and GTR + G + invari- 
able sites (GTR +I+ G) for dataset 2, as selected by 
jJModelTest [40], and both runs under a lognormal relaxed 
clock and constant population size. For dataset 1, the 
MCMC chain ran for 20,000,000 steps sampling every 2000 
steps and for dataset 2 for 100,000,000 steps sampling every 
5000 steps. The phylogenetic trees were summarized using 
Tree Annotator 1.8 and edited using iTOL 4.2.3 [41-43]. 

Association between host traits and tree topology was sta- 
tistically tested based on the posterior set of trees (PST) gen- 
erated by BEAST, using the software BaTS 0.9 [44]. The 
Parsimony Score (PS), Association Index (AI), and 
Monophyletic Clade (MC) metrics were analyzed, and values 
of p< 0.01 were considered significant. 
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Images Design and Attribution 


The map presented in Fig. 1 was designed by the authors using 
QGIS 3.2.3 based on maps of cities and states provided by 
Brazilian Institute of Geography and Statistics (IBGE) [45, 
46]. Vectorial images of animals used in Figs. 1 and 2 were 
provided by vecteezy.com, www.freepik.es, and pixabay.com. 


Results 
Samples and Coronavirus Detection 


A total of 103 bats were collected, 64 from Sao José do Rio 
Preto region and 39 from Barreiras region. In both locations, 
males were more sampled than females, 60.9% and 66.7% 
respectively. In Sao José do Rio Preto, the species that was 
most sampled was Molossus rufus (21.9%) followed by 
Molossus molossus and Artibeus literatus (20.3% each). In 
Barreiras, Molossus molossus (28.2%) was most sampled 
followed by Artibeus literatus and Artibeus planirostris 
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Fig. 1 Map of Brazil showing the locations where infected bats where 


collected. The host families of infected animals are represented in the 
respective sampling location (see legend). São Paulo state—dark blue; 


São José do Rio Preto 
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(23.1% each). Online Resource 3 present the information on 
sex, family, species, and organs of all animals collected. 

The presence of CoV RNA was tested by reverse transcrip- 
tion NESTED-PCR specific for the nsP12 genomic region. 
From bats collected in São José do Rio Preto region intestines, 
liver, and lungs were tested. As for the animals from Barreiras 
region, only intestines and liver were tested, due to availabil- 
ity. Positive samples were sequenced and the resulting se- 
quences were submitted to BLASTn (https://www.ncbi.nlm. 
nih.gov/BLAST/) in order to confirm CoV identification. The 
results revealed that 12 animals were infected with 
coronavirus, 11 from Sao José do Rio Preto and 1 from 
Barreiras. Viral RNA was detected in more than one organ 
in some animals totalizing 17 positive samples. Only one 
animal, SJRP03, tested positive in all organs. Coronavirus 
was mostly present in intestine samples, except for SJRP06 
and SJRP10. Table 1 displays information on family, species, 
gender, and the organs where the virus was detected. 

Considering the regions separately, the prevalence of CoV 
in bats from Sao José do Rio Preto region was 17.18% (11/64) 
and 2.56% (1/39) in bats from Barreiras region. Also, males 
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Sao José do Rio Preto region—light blue; Bahia state—orange; Barreiras 
region—light orange 
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Fig. 2 Bayesian phylogenetic tree based on a dataset of 52 nucleotide 
sequences of 426 bp from the nsP12 genomic region of CoV. Nodes 
display posterior probabilities higher than 0.7. The analysis includes 17 
sequences from this study along with representatives of the four CoV 


were more infected than females since in São José do Rio 
Preto, CoV was detected in 23.1% (9/39) of males and in only 
8% (2/25) of females. In Barreiras—BA, the only bat infected 
was male. The locations where infected bats were collected 
are presented in Fig. 1. It also show the families of bats that 
hosted CoV found in each location. 


Phylogenetic Analysis 


A total of 17 sequences, from 12 different animals, were gener- 
ated in this study. The first analysis we performed was to deter- 
mine the phylogenetic relationship of the sequences identified in 
both regions with representatives of the four Coronavirus genera. 
Based on dataset 1, composed by 52 nucleotide sequences of 
426 bp (see “Methods”), we ran a Bayesian phylogenetic anal- 
ysis. The phylogenetic tree indicates that the 17 sequences from 
this study belong to the -CoV genus since they group in a 
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genera, represented by different colors according to the legend. Images 
beside the tips of the tree are a graphic representation of the host. 
Accession numbers of each sequence is presented in Online Resource 
1. Sequences generated in this study 


monophyletic branch together with all sequences from this ge- 
nus, with strong branch support (1.00) (Fig. 2). 

We decided to analyze our sequences in the context of Bat 
Coronaviruses from the a-CoV genus. To this end, we assembled 
dataset 2 composed by 139 nucleotide sequences of 338 bp of the 
nsP12 region, which included all sequences from this group 
available in GenBank, plus our sequences. A Bayesian MCMC 
analysis using BEAST 1.8 was performed. The summarized tree 
is presented in Fig. 3. The analysis of the phylogenetic tree re- 
vealed a tendency of viruses isolated from the same host family 
to group in monophyletic branches with strong branch support, 
although there were some exceptions. In order to test for this 
association, we performed a statistical analysis, based on the 
posterior set of trees using BaTS software. According to the 
results, there is a positive association between the clustering in 
the phylogenetic tree and bat families with a p value < 0.01 for 
PS and AI. As for MC values, all families displayed significant 
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Table 1 Specimen, family, species, gender, and the organs that tested positive for Coronavirus 

Bat ID Family Species Gender Intestines Liver Lungs 

BtCoV_ Molossidae Molossus rufus Male + + + 
SJRP03 

BtCoV_ Molossidae Molossus rufus Male + = = 
SJRP04 

BtCoV_ Molossidae Molossus rufus Male + + - 
SJRPO5S 

BtCoV_ Vespertilionidae Eptesicus sp Female — + + 
SJRP06 

BtCoV_ Phyllostomidae Phyllostomus discolor Male + = 7 
SJRPO8 

BtCoV_ Phyllostomidae Glossophaga soricina Male z = + 
SJRP10 

BtCoV_ Molossidae Molossus molossus Female + - -= 
SJRP16 

BtCoV_ Phyllostomidae Phyllostomus discolor Male + — — 
SJRP22 

BtCoV_ Phyllostomidae Artibeus literatus Male + S = 
SJRP44 

BtCoV_ Molossidae Molossus rufus Male + E T 
SJRP49 

BtCoV_ Molossidae Molossus rufus Male + + = 
SJRP52 

BtCoV_BA191I Molossidae Molossus molossus Male + -— - 


association except for Mormoopidae, since there is only one 
virus identified in this family in our analysis (Table 2) (Online 
Resource 4A). 

We also tested for an association between phylogeny and 
the location where the virus was sampled. The results revealed 
a positive association between this trait and phylogenetic clus- 
tering based on AI and PS metrics. The MC metrics was sig- 
nificant for all locations with more than two representatives 
(Table 3) (Online Resource 4B). (Sampling site of each 
sequence used in the analysis is provided in Online 
Resource 2). 


Discussion 


The outbreaks of the G-CoV SARS-CoV and MERS-CoV 
evidenced the importance of bats as natural reservoir and 
as a source of spillover of Coronaviruses into the human 
population [26]. Understanding the diversity and the eco- 
logical aspects of these viruses in bats is important to pre- 
vent new outbreaks. Here, we investigated the presence of 
Coronavirus RNA in bats, living in urban and peri-urban 
areas, from two geographically distant regions in Brazil, 
Sao José do Rio Preto—SP and Barreiras—BA. Our results 
revealed the presence of Bat Coronavirus from the ~-CoV 
genus in specimens from both regions. None of them are 
closely related to pathogenic viruses to humans. A preva- 
lence of 17.2% was observed among bats sampled from the 
Sao José do Rio Preto region which is similar to the rate of 


16.7% reported in bats from Hong Kong [47]. Although 
one study in the UK reported a comparably high BtCoV 
prevalence of 23% [48], reported prevalence rates usually 
range from 3% to 10% [13, 16, 17, 22, 24, 49-53], similar 
to our current findings in the Barreiras region (2.56%). A 
global study with over 12,000 samples reported a BtCoV 
prevalence of 8.6% [54]. The CoV prevalence among bats 
from the Sao José do Rio Preto region is therefore higher 
than estimated for Barreiras as well as compared to other 
studies. Several factors might influence the estimates, like 
the number of specimens sampled, host species diversity, 
number of sites of collection, geographic regions, among 
others [54]. In this study, 54.5% (6/11) of the infected 
specimens from Sao José do Rio Preto region were collect- 
ed in the same site. The location is a private property 
known to harbor about six different species of bats in the 
roof of the house. The fact that these animals were collect- 
ed in the same location might have influenced the high 
CoV prevalence in this region, since they live in close 
contact aiding the transmission. Another factor that might 
have influenced is that in addition to intestines, we also 
analyzed lungs and liver for the presence of CoV, while 
most prevalence studies only test fecal or intestine sam- 
ples. Although most of the infected animals tested positive 
in the intestines (81.8%—9/11), in two animals, CoV was 
only detected in other organs (lungs and liver), raising the 
chances of viral detection. Other authors have also reported 
the presence of CoV RNA in tissues and secretions other 
than intestines or fecal samples, such as kidney, heart, 
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Fig. 3 Bayesian phylogenetic tree based on a dataset of 139 nucleotide 
sequences of 338 bp from the nsP12 genomic region of CoV . The 
analysis includes sequences of BtCoV from the a-CoV genus. Nodes 


lungs, and oropharyngeal swabs [25, 55, 56]. However, 
this is the first report of BtCoV detection in the liver. 
Since the animals presented no clinical symptoms, the sig- 
nificance of the detection of virus in lungs and liver re- 
mains unclear. Experimental infection of bats with CoV 
from the Betacoronavirus genus has been tried by 


Table 2 Statistical 


analysis of the Statistic Bat family p value 
association of host 
family and phylogenetic Al 0.00 
clustering performed on PS 0.00 
Pa MC Miniopteridae 0.01 
Parsimony Score (PS), f 
Association Index (AT), MC Molossidae 0.01 
Monophyletic Clade MC Vespertilionidae 0.01 
(MC). Values ofp<0.01 MC Phyllostomidae 0.01 
ge considered MC Rhinolophidae 0.01 
significant 
MC Mormoopidae 1.0 
MC Hipposideridae 0.01 
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KC110785 
KC110779 


Legend 

Molossidae 
Vespertilionidae 
Phyllostomidae 
Hipposideridae 
Miniopteridae 


Rhinolophidae 


Mormoopidae 


display posterior probabilities higher than 0.7. Tips are colored 
according to the host family (see legend). Sequences generated in this 
study 


Watanabe and colleagues [57]. Results revealed that the 
virus was only detected in the intestines. In our study, the 
infection of lungs and liver was not restricted to one host 
species or family since it was found in Molossus rufus 
(Molossidae), Eptesicus sp. (Vespertilionidae), and 
Glossophaga soricina (Phyllostomidae). 

We analyzed the phylogenetic relationship of the virus 
identified in this study with other BtCoV from the -CoV 
genus (Fig. 3). The sequences generated in our study are 
distributed in three distinct monophyletic branches. Most 
sequences, 13 out of 16, of CoV identified in bats from Sao 
José do Rio Preto grouped into a monophyletic branch with 
strong posterior probability (0.77). In the same branch, 
there are also three sequences of CoV identified in bats 
from the cities of Araçatuba (KC886321) and Penapolis 
(KT717388, KT717389), about 100 km distant from Sao 
José do Rio Preto. This cluster shares a most recent com- 
mon ancestor, also with high branch support, with the virus 
we identified in a Molossus molossus bat in Barreiras 
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Table 3 Statistical analysis of the association of location and 
phylogenetic clustering performed on BaTS software. Parsimony Score 
(PS), Association Index (AI), Monophyletic Clade (MC). Values of p < 
0.01 were considered significant 


Statistic Sampling location p value 
Al 0.00 
PS 0.00 
MC Hong_Kong 1.00 
MC Barreiras BA 1.0 
MC S.J.Rio Preto_SP 0.01 
MC Bulgaria 0.01 
MC Panama 1.0 
MC Porto Alegre RS 0.01 
MC CostaRica_ 1.0 
Tirimbina 

MC CostaRica_SanJose 1.0 
MC Aragatuba 1.0 
MC Mexico 1.0 
MC China_Henan 1.0 
MC CostaRica_Parrita 1.00 
MC Mojiang 0.01 
MC Ghana 0.01 
MC Fénix_PR 1.0 
MC Iguagu_PR 1.00 
MC Ivai_PR 0.01 
MC Andradina_SP 1.0 
MC Penápolis _SP 0.01 
MC SudMennucci_SP 1.0 
MC Campinas_SP 0.01 
MC Ribeirao_Preto_SP 1.0 
MC Valinhos_SP 1.0 
MC Sumaré_SP 1.0 
MC China_Guizhou 0.01 
MC China_Guangdong 0.01 
MC Hainan_China 1.0 


(BtCoV_BA19I). The sequence BtCoV_SJRP44I groups 
apart from all other in a monophyletic branch (0.99) to- 
gether with sequence KT717382. Both viruses were detect- 
ed in Artibeus literatus bats. It is interesting that 
KT717382 was identified in Iguacu National Park—PR, 
about 850 km from São José do Rio Preto—SP, in 2012 
[22]. Sequences BtCoV_SJRPO8I and BtCoV_SJRP22I, 
both from Phyllostomus discolor bats, group together in 
another monophyletic branch sharing a most recent com- 
mon ancestor with a CoV from a P. discolor (JQ731783) 
bat collected in Panama [23]. Both clusters include all 
strains, included in the analysis, identified in the respective 
species. This indicates that the BtCoV identified in these 
hosts could be species-specific, being more closely related 
to the viruses identified in the same host, even if they are 


geographically distant. In fact, other sequences from our 
study, from hosts collected in the same sites, were infected 
with BtCoV that grouped apart from the P discolor and 
A. literatus clusters (Fig. 3; Online Resource 2). An alter- 
native explanation is that the occurrence of related strains 
in distant locations is a consequence of transmission during 
bat migration. Given the genetic distance between the 
strains sampled, this would obviously not be due to direct 
transmission but would presumably involve unsampled in- 
termediate strains. Several bat species are known to mi- 
grate seasonally following resources, such as food. 
Migration studies report that small- and medium-sized bats 
can migrate long distances, with a maximum migration 
distance of 1905 km reported for a Vespertilionidae bat 
[58]. However, there are no studies of migration patterns 
for P. discolor and few for A. literatus, none of which 
report such long distance [59, 60]. It is therefore more 
likely that the close relatedness of BtCoV from 
P. discolor and from A. literatus is a reflection of the virus 
“host specificity.” 

Overall, the BtCoV phylogenetic tree shows highly sup- 
ported branches clustering primarily according to host 
family. Some families have two different clusters like 
Molossidae, Phyllostomidae, and Rhinolophidae while 
Miniopteridae and Hipposideridae present only one each. 
Since there are some exceptions and CoV from 
Vespertilionidae bats does not seem to group together, we 
statistically tested the hypothesis of a positive relationship 
between host family and phylogeny and it was confirmed. 
However, due to the exceptions, we decided to test for a 
possible association with location where the hosts were 
collected. For example, as previously described, most 
CoV from Sao José do Rio Preto grouped into a monophy- 
letic and most of them belong to the Molossidae family. 
Still, among them there are sequences from virus identified 
in two animals from other families, Vespertilionidae 
(BtCoV_SJRP06) and Phyllostomidae (BtCoV_SJRP10). 
So, we statistically tested for an association between the 
location where the host was collected and phylogeny and 
there is also a positive association. Our results indicate that 
despite the association between BtCoV and host family 
seen here, the close contact between animals from different 
families may result in virus transmission. We conclude that 
both traits, host family and location, influence BtCoV phy- 
logenetic relationship. 

Here, we report the detection of Alphacoronavirus in bats 
from Brazil, not only in intestines but also in lungs and liver. 
We found a higher prevalence of BtCoV in animals collected 
in Sao José do Rio Preto—SP region, when compared to 
Barreiras—BA. The investigation of CoV in bats throughout 
the world remains important to characterize and understand 
the circulation of these viruses, identifying potential sources 
of spillover of pathogenic strains into the human population. 
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